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Abstract 

A spherically symmetric charged ideal fluid solution of Einstein field 
equation is given in the presence of the cosmological constant and 
two well known example of this type of solution is presented. If the 
matter is confined in a region, the exterior spacetime is considered as 
RN-de Sitter (Reissner-Nordstrom de Sitter) and to complete solution 
matching conditions are examined. We show that the function which 
is related to the dynamics of the system will determine the fate of 
the system: expansion, contraction or bouncing situations may oc- 
cur for different configurations. The initial conditions of the matter 
determine the final form of the system and therefore the nature of 
the singularities in the presence of the electric charge and the cosmo- 
logical constant is examined to reveal their effects on the singularity 
formation during collapse. 

PACS number(s): 04.70.Bw, 04.20.Dw, 04.20.-q 
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1 Introduction 



Cosmological observations indicate that the expansion of the universe is ac- 
celerating. A simple phenomenological interpretation of the data in terms of 
dark energy has been successful thus far [1]; dark energy is assumed to have 
negative pressure, unlike ordinary mass-energy, and thus leads to a negative 
force that may account for the acceleration of the universe. Recent obser- 
vations indicate that universe is dominated by dark energy (~ 70%), which 
can be thought of as a perfect fluid with an energy-momentum tensor given 



where Hd+Pd — and < 0. The inclusion of source term (1) in Einstein's 
field equations amounts simply to a cosmological constant A given by 
A = —pr) > 0. In this way, dark energy may be represented by a positive 
cosmological constant A in the Einstein field equations 



The cosmological constant has a long history in Newtonian and relativistic 
gravity theories. It was introduced into general relativity by Einstein as a 
means of balancing the gravitational attraction of the matter on cosmological 
scales, leading to the Einstein static universe model. Alternatively, it can be 
thought of as a measure of the energy density of the vacuum. 

In this paper the gravitational motion of charged matter is considered 
in the presence of a cosmological constant. To simplify matters, a perfect 
fluid distribution is considered that is electrically charged and undergoes 
spherically symmetric collapse or expansion in the absence of external fields. 



by 



T^''^{i^d+Pd)u''u^+Pd9'"', 



(1) 




(2) 
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The exterior field is given by the Reissner-Nordstrom-de Sitter metric 



where M is the net mass and Q is the net charge of the system and A is the 
cosmological constant of the RN-de Sitter spacetime. 

The gravitational collapse phenomena is still open problem in the general 
relativity and it has not been taken its final form yet. Throughout the studies 
it is shown that the initial conditions arc determinative of the final fate of 
the collapse [2]. If the collapse can not be stopped in an equilibrium state 
and allows to formation of the singularities as a result, the end state of 
the collapse can either be black hole or naked singularity depending on the 
character of the singularities. If all singularities are hidden behind an event 
horizon singularities can not be seen by distant observer and the end state 
becomes " a black hole" , or if the singularities are bare and can be visible by 
distant observer, the end state becomes "naked singularity" [2]. 

One of the physical features affected on the gravitational collapse is shear. 
In [3], shear effects on the gravitational collapse of the spherical massive cloud 
with non-radial pressure are studied and shown that sufficiently strong shear 
effects near singularity delay the formation of the apparent horizon and allow 
the formation of the naked singularity. 

In [4] effects of the cosmological constant on the gravitational collapse 
of the pressureless matter is studied and shown that positive gravitational 
constant plays repulsive role and slows down the collapse process. 

As a property of matter how does electric charge affect the collapse phe- 
nomenon? This question finds some answers through the following papers 
and needs to be completed. A solution for the spherical symmetric charged 



ds^ = -(1 : — \—:^~Trr )dt 

r r 6 



(1 




^f''{d0''+sm''§d(f)^) 



(3) 
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stars are considered in point of view formation of the black holes and voids 
in [5]. In [6] , the gravitational collapse of the spherical symmetric charged 
radiating Vaidya-RN type spacetimes are studied. To avoid singularities due 
to charge of the matter in spherical symmetric collapse is examined in [7] . As 
a stellar model, relativistic structure, stability and gravitational collapse of 
the charged fluid is studied in [8] , the specific values of the electric charge of 
the fluid allows formation of naked singularity besides black hole formation. 
A spherically symmetric charged ideal fluid is examined in [9] due course of 
gravitational collapse. Similar to solution presented in [9], in this work, we 
give a solution of Einstein field equations in dc Sitter spacetime which is in 
isotropic form. Furthermore, we examine the effects of the electric charge on 
the collapse phenomenon in the presence of the cosmological constant. 

In the following section, a spherically symmetric solution of the Einstein- 
Maxwell equations is given in the presence of the cosmological constant. The 
exterior spacetime of the charged fluid sphere is considered the RN-de Sitter 
spacetime and matching conditions about two distinct (interior and exterior) 
regions are examined in section 3. Two well known examples of this type 
of solutions, RN-dc Sitter and Mc Vittic-dc Sitter arc given. Section 4 is 
devoted to the gravitational collapse of the charged fluid and formation of 
the singularities. Since the co-moving character of the spacetime may give 
coordinate dependent results, it is necessary to use coordinate free "null 
geodesies method" where the nature of the singularities does not change its 
character. Therefore, the nature of the singularities are investigated null 
geodesic method. Throughout the paper we use units such that c — 1 and 
(SttG^ 1). 
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2 Interior solution 



Imagine a co- moving system of coordinate for the interior {t, p,9,(f)) that 
remains at rest with the moving charged matter and is given by 

ds'' = -a^ + dp'' + R' p'dO^ + p' sm' 0d<l>'' (4) 

where a — a{t, p), b — b{t, p) and R — R{t, p) are arbitrary positive functions 
of time coordinate t and radial coordinate p. To consider shear free motion 
of the matter that is the ideal fluid, we take R — ph{t,p). (4) is assumed 
to be the solution of the field equations (2) with a cosmological constant Aq 
and a source T^^ = + T^Jf , where 

+ P) u^u^ + p g„u (5) 
and the electromagnetic energy-momentum tensor is defined by 

= Sl^-'^F^^F^. - \g,.F^pF"(^) . (6) 

Here electromagnetic field tensor is F^i, — d^^A^ — d^A^ and A^ is the vector 
potential. The spherical symmetry of the spacetime ensures existence of a 
radial electric field in general. 

After a suitable choice of the gauge is Ai = 0, At = ^{t,p) chosen, ac- 
cording to the metric (4) the only non-zero component of F^^ becomes 
Ftp — —d^/dp. Then, non- vanishing components of the electromagnetic 
energy momentum tensor for the spacetime with the gauge chosen above are 
given by 

* pH^' 
Tl = T^P^-T,' = -T} (7) 
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where q is the total charge of the fluid. It is easily seen that T = trT^ = Q 
as expected. The Maxwell equations 

d , 



gF''^] = An^r (8) 



with — ab^p^ sin 9 for the spacetime (4) become 

It is seen that total charge q defined by 

6 9$ q 



2 (10) 

a op 



is independent of time and related to the charge density ( by the equation 

Wb\. (11) 



dp 

The electric current is given by J* = C'^* ™ terms of the electric charge density 
and 4- velocity — a~^. Since the total energy momentum tensor is the sum 
of the electromagnetic T^Jf and matter part T^, the total energy- momentum 
components will be obtained as 



where p are p are the matter-energy density and pressure respectively. It is 
supposed that the matter field satisfies the weak energy condition. For any 
timelike vector v'^ 

T^^v^v^ > (13) 

which gives 

p > 0, p + p > 0. (14) 
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Furthermore, the conservation of the energy u^T'^''.^ = gives the relation 

1 dfj, _ 3 db 
(l^ + p) dt ~ b dt ' ^ ' 

Now we are going to examine the geometry part of the problem given 

by the metric (4) and consider the Einstein tensor G^i, — R^i, — Rg^i,/2 as 

follows 



Gtp — 



62 ^2 



-2a X 



b \blpb 



ab , 



G 

Ge 

Ga, 



pp 



+ 2 



a'b' 2 (y a'' 



+ 



+ 



a b p \b a 



^ f 2 ^ + - - 2 ^ 
a? \ b a b , 



1 /a' b'^ 
p\a b ^ 



62 



ab 



' 2": 

b. 



sin^ e Ga 



(16) 



where "dot" and "prime" represent derivatives with respect to t and p, re- 
spectively. 

By using the energy-momentum tensor of the fluid obtained in (12), the 
field equations (2) take the form 



Gtp = 0, 



Gu - An 



— Gpp+Ao = p 
Gpp 5- Gee — 



pH^ ' 
2g2 



(17) 
(18) 

(19) 

(20) 



where Aq is the cosmological constant for the interior region. For the sake of 
the generahty, the cosmological constant of the interior region Aq is taken to 
be different from the cosmological constant of the exterior region A in the 
beginning. 

Then, the trace of field equations becomes 

_ 2 /a" ^ a'b' 6'2 ^ 2b" ^ 2a' ^ 4b'\ 6 /ab 6^ ^\ 
b"^ \ a ab 6^ b ap bp J a"^ \ab b'^ bj 



Equation (17) has the solution 

b^abk{t), (22) 
where k{t) is an arbitrary function of time, then equation (20) becomes 



'a" b"\ (I ^b'\(a' b'\ 2g2 



A solution of the full Einstein field equations can be found by following the 
method given in the reference [9] as follows: 

■r2 , 1 An A I/A^, 

1 + - + ^ 



where^=^(^l^^j, A(0 = , W ^ {a - ^r'){5r' - 

and {a5 — (5^) > 0. a, /3, 5, 7, 770 are all real and Xq > 0, u > are pos- 
itive constants. Moreover, f{t) is positive arbitrary function of time and 

cosmological constant Aq . 
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Since the radial coordinate transformation does not change the co-moving 
character of the metric, for the sake of the brevity, we used the transformation 

1/2 

/ ucp" -|- jj \ 

r 



Then, the physical quantities mass-energy density and pressure satisfying 
relations (18), (19), can be written in following form 

// = -Ao + 3(4) + ^ -3 L 



p = Ao 



J) K 

r'iHvi - A^) - 8rAo/ + Ar'f) f 2 iVp / 
PN, fN^ 



with 

No = -vo' + Ao' + 4rAo/ + 4rV' , N, = rjl - Xl + Ar'f , 

r^^yp'j — ao) ^ 
+{f3-f-a 5)V {-aPvo - {'2 f r + X^f {A f ^ 5 - al3 Aq') 
+770^ (-4/2 P{-Za^ + P l3^ + PaS) + Sfra^Xo + ^apXo^) 
+77o'(2/r + Ao) (8/V^(-/37 - a5 + 3 7 5) 

+4 f r\-4ap + Pf^'y + PaS + r V) Ao - 10 / r a /? Aq^ - 3 a /3 Aq^' 

(27) 

and the total charge of the fluid is given by 



^0 



[a — 7r 

)3/2(5^2_^)3/2_ (28) 



(26) 



Consider the spacetime given by the hne element (24) describes whole space- 
time then, by setting of arbitrary constants {a. j3. 5, 7, rjo, Aq, v) and func- 
tion of time / in special forms, well known solutions of Einstein field equations 
such as RN-de Sitter and charged McVittie solutions can be obtained. 



1) RN-de sitter solution. 

The line element of the isotropic RN-de Sitter spacetime in [10] is given by 



1 - 



m 



+ 



1 + — ^ 

wr 



1 + — ^ 

wr 



yj2j.2 



{dr' + r'de'' + r' Sin'' e'- 



(29) 

where m is the mass and q is the electric charge of the black hole and w — 
is function of the Hubble parameter which is in general a function of time. 

If we take 

f^w^ *, Ao = 2m, rjo ^ 2q, a = 5 = 1, /3 = 7 = (30) 

the line element (24) reduces to isotropic RN-de Sitter line clement (29). 

In this configuration rjo/2 and Ao/2 are considered as total charge and 
mass of the black hole, respectively. Since f / f — ^Ao/3 > , therefore 
b/b > 0, then the spacetime expands in time and H — const. — y^Ao/S 
corresponds to the cosmological constant of RN-de Sitter spacetime. Fur- 
thermore, since m (or Aq) and q (or 770) are constants, /i — p — 0, are all zero 
as expected. 

The spacetimc (24) contracts only if /// is negative. If we take / = e~^^°^^ *, 
it will correspond to isotropic RN black hole in anti-de Sitter spacetime. Here 
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the physical quantities mass-energy density and the pressure are also zero as 
the previous case. 

The cosmological constant is the reason of the time evolution of the prob- 
lem. Therefore, if we take Aq = 0, the contraction (expansion) of the space- 
time disappears, it becomes static. 

2) Chcirged McVittie-de Sitter solution. 

A perfect fluid solution of Einsteins equations corresponding to the Schwarzschild 
field embedded in a RobertsonWalker background is given by McVittie [11]. 
In expanding universe the McVittie solution represents a white hole, whereas 
it represents a black hole in contracting universe [12]. In addition to the re- 
pulsive effect of the expansion of the spacetime, the repulsive character of the 
electric charges (Coulomb force) of the fluid will be greater in small regions 
and will support formation of white hole in the beginning in the McVittie-de 
Sitter spacetime. 

If we consider (r/o 7^ /? = 0, a and S positive) in (24), the charged 
Mc Vittie-de Sitter solution is obtained and the isotropic metric components 
become 

iar'f + {S + ^r') {rjo' - Xo') 



a — 



Aar^P + 4r ^a{5 + 7 r^) / Ao - (5 + 7 r') (r^o" - Ao') ' 



^ ^ 4 r + 4 r ^a(^ + 7 r^) / Ao - (^ + 7 r^) ivo' - Ao') 
The matter density and the pressure can be written as 

128/3 (2/(0; - r^^)%o'^ + 3r^a'^j5^{2fr + Xof + r^a^So" {V (« - ^r'^l) - 3r7Ao)) 

aM2(V-(2/r + Ao)')' 
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P = Ao 



256 f^d-fr 





and the charge becomes 

(33) 

By taking cosmological constant Aq = 0, and electric charge 770 = 0, un- 
charged ordinary McVittie solution can be recovered 

Mass-energy density and the pressure of the fluid become 




[V 384Ap^7r5/3 

j) + (2r/ - Ap)3(2r/ + Ap)3 ' 



256 (5 7 /A 5Ao + 8Aor/-4r^ f _ /A (2r/ + 

(2r/ - Ao)(2r/ + Ap)^ + V/^J (2r/ - Ao)(2r/ + Ap) l^/J "(^7/^ 



Now let us consider the charged fluid is confined in a region in which dynam- 
ical evolution of the system is described by the time dependent function f{t) 
and the exterior spacetime to the confined matter is RN-de Sitter. Exterior 
and interior regions separate the spacetime into two distinct parts such that 
they meet on the boundary surface. To complete the solution, it is neces- 
sary to show that the distinct solutions must satisfy boundary conditions on 
the boundary surface. In the following section we will investigate matching 
conditions. 
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3 Matching conditions 



Let us consider a spherical boundary surface which divides spacetime into 
two distinct four-dimensional manifolds which admit S as their boundaries 
at rs = = "Te = const.. 

Let t, f, 6, (j) be the Reissner-Nordstrom de Sitter coordinates for the 
matter-free region. Then, the metric is 

where 

r r'^ 6 

with the interior metric is in original form and given by 

ds'^ = -a^ de + dp" + 6^ p2 (c^^2 ^ g.^2 Q ^^2^ (37^) 

where a and h are functions of coordinates p and t. Under the coordinate 
t r ansf ormations 

e^e, = (38) 

motion of the boundary surface can be given by the following equations 

/+:r-fs(t) = 0, r:p-ps = 0. (39) 

To match the exterior spacetime with the interior spacetime we use the Israel 
junction conditions [13]. These conditions require the interior and exterior 
solutions of the gravitational field equations to be joined smoothly up to a 
coordinate transformation but the partial derivatives may change discontin- 
uously across the boundary surface of the matter. 
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Let (is| be the line element of the boundary surface S, ds\ represent the 
exterior and ds'^_ represent the interior spacetimc line elements. The junction 
conditions which state the equality of the first fundamental forms and the 
discontinuity of the second fundamental forms can be given as 

ds| = o?s!_ Is = rfs+ls , 

Kt,-K-^- Qij K = T,,. (40) 

where K — g'^^ {Kij — K^^ and Tij is the surface energy momentum tensor. 
In case of vanishing surface energy momentum tensor i.e. for = 0, the 
discontinuity condition reduces to the equality of forms i.e., equality of the 
extrinsic curvatures 

The continuity of the metric components in (40) at p = ps on the boundary 
surface gives the following relations 

= Ps&(t,Ps) , dr = a{t,pj:)dt, (42) 

dr 



t^'^'>-Ak) '''' 



and 



df db p db . 
dr dr a dt ^ 



'drV 



\2 



(45) 



\dt J A + {dr/dT] 

The non-zero second fundamental forms for the interior and the exterior 
regions are given by 

K^, = smHK^^ = p{bpy (46) 
14 



1 ( dr^ 



(47) 



and 



- (48) 
a 

d^d^~ d^d^^2A~df [d^J d^~2W\d^) ' 
By using the equality of the angular components of the second fundamental 
forms (46, 47) we obtain the following relations 

Ar^^pibp)', (50) 

Since b is function of time dependent function /, (51) gives a condition that 
/ should satisfy at any time t . Furthermore, the equality of the timelike 
components of the extrinsic curvatures, (48) with (49), gives 

{Adi/dT)^r 1 da . . 

df/dr ah dp ' 

in other words it corresponds to Gtp = that is, no new information is 
obtained from timelike components of K'^^^. Since the radial pressure of the 
fluid is zero on the boundary surface, this condition reduces to the continuity 
of the energy momentum tensors in radial direction — In which gives 

p{t,p)-Ao^^^^^-^-A. (53) 

For the sake of the generality we started with different cosmological constants 
for interior and exterior regions, but the continuity of the energy momentum 
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tensor in the radial direction compels their equality "Aq = A" . Then, the 
electric charge distribution can be written as 

where Q is the total electric charge of the fluid confined in the region p < ps, 
and Wj^ is the value of at p = ps given by (24). 

4 Gravitational collapse 

If the collapse phenomena allows their formation two types of singularities 
may form during collapse: physical and spacetime singularities [2]. Physical 
singularities make physical quantities (such as mass-energy density, pressure) 
singular and the space-time singularities make the metric components and 
the curvature indefinite. In the gravitational collapse manner, among the 
spacetime singularities the shell focusing and shell crossing singularities are 
being considered. The shell crossing singularity occurs at distances where 
change of the radius of the fluid sphere in radial direction is zero R' = 
(with R > 0), and the shell focusing singularity forms at distances which 
make radius of the fluid sphere zero (R ^ 0). The shell crossing singularities 
can be considered weak with respect to the shell focusing singularity in the 
gravitational collapse treatment [14]. Therefore, we are only interested in 
the formation of the shell focusing singularity, i.e., i? — > as p — > with 
R'>0. 

In the literature many factors effective on the formation of the naked 
singularities are examined [2], [3], [15] and it is pointed out that, one way 
of a singularity to be naked is to disturb the apparent horizon surface and 
delay its formation [3]. According to the idea, if the trapped surface forms 
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before the singularity surface then, the singularity becomes hidden inside a 
black hole. Otherwise, the trapped surface forms after the singularity surface 
and the singularity becomes naked. In another words, if the time period for 
the formation of the event horizon is longer than the time period for the 
formation of the singularities, singularities become bare and they can be 
seen by distant observer. This criterion is probably easy and efficient, but 
it is not clear if this is always equivalent to naked singularity formation and 
in a way it is coordinate dependent statement. In one coordinate system 
these two timings may be related in a certain way, but may not be related in 
another coordinate system. Therefore, the coordinate independent and a full 
proof condition "the families of null geodesies come out of the singularity" 
should be examined [2]. We will give this analysis in the next section. 

Let us write metric components of (24) in the radial coordinates p ex- 
plicitly to examine under which circumstances the situation corresponds to 
collapse, expansion or bounce 

Ob 



4/2( V + /5) + 4/ AoV( V + + ^) + (Ag - nl) ilp' + 8) 

_ Apjap^ + /3) + 4/ Ao7(ap^ + (3){^(^ + 5) + (Ag - 4) (tp^ + 5) 

(54) 

and define the physical radius R, the radius of 2-sphere 



Afjap' + /3) + 4/ AoV^( V + (3)i^p^ + 6) + (Ag - yj) j^p' + S) 
^ ^ Af{ap^ + f3){^p' + d) 

(55) 

As the radial coordinate p — > the physical radius shrinks to zero. We know 
that if the spacetime is in the isotropic form the shear tensor is zero then it 
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reduces to 

^^ (-) 

By means of shear free property of the spacetime, the change of rate of the 
physical radius that is, the expansion rate of the spacetime is given by 

This equahty states the dynamics of this isotropic collapse problem is related 
to ratio ///, and collapse (expansion) situation can only occur for the neg- 
ative (positive) values of ///. Since / is a positive function of time, only if 
/ < (> 0) the negative (positive) expansion rate corresponds to contract- 
ing (expanding) physical radius R . If / < changes its sign after a time 
period then, it is called bounce. / = will be static solution. The time 
evolution of the physical radius is explicitly written by 

R^oJ^t 4/ W + /^) - (Ag - ^g) (7^^ + ^^) 

R f f AP {ap^ + /3) + 4/ Ao^( V + p)(^p2 + S) + (Ag - rj^) {^p^ + 5) ' 

(58) 

from equations (56) and (57). By analyzing the time dependency of the R we 
see that four different situation can be obtained: collapse, expansion, stable 
and bouncing cases. 

In the literature it is shown that positive cosmological constant delays 
the formation of the singularities [16, 4]. By considering slowing down effect 
of the positive cosmological constant on the collapse process let us take time 
dependent function as f — e~^^~^^°^^^* where c is a positive constant. Since 
it/R = -(c - 0V3) a, with -y/g^ — a > 0, the expansion (contraction) 
becomes dependent directly to the values of the cosmological constant Aq. If 
3c^ > Aq the collapse will continue to a certain radius and be stopped by the 
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Coulomb's repulsive force by the charge and it will reach to the central point 
for uncharged matter. If 3c^ < Aq radius will expand until it reaches to the 
boundary of the exterior region that is to the apparent horizon stated by the 
matching condition (51) . Purtheremore 3c^ = Aq corresponds to static case. 
In this example all three possible situations are obtained: expansion, crunch 
and stable cases. 

If the cosmological constant effects are dominant in the dynamics, it will 
be more convenient to choose sample function as / = e"^*~'\/^°/^*^ > 0) to 
emphasize the cosmological constant dependency. In this model the change 
of the radius with time R/R = (— c + 2Aot)a, a > 0, is positive in the 
beginning and negative for late time t. It means that the radius of the fluid 
will decrease with time and after a period of time, here for t > c^Js/AAq , the 
radius will start to increase so, "bounce" situation is obtained, a — case 
will be examined in the following part in details. 

Physical singuleirities. 

When pressures are non-zero, dynamical evolutions, as allowed by the Ein- 
stein equations, are equally important as the initial data is to determine the 
final fate of collapse. [17] . Dust solution of the gravitational problem is 
highly important, but the isotropic form of the spacetime and set up of the 
problem does not give the dust solution consistent with the conditions. 

When their singularities are examined, it is seen that matter and charge 
densities (26, 28) arc regular everywhere but the pressure which is subjected 
to the the weak energy condition (14) and satisfying the energy conservation 
condition (15) 

P^-if^+l^^) (59) 
3a/ 

diverges at the distance p = Ps which makes the metric component 
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' -4/3/2 + 5(Ag-77o^) 



4/^(V + /3)-(A^-r/o^)(7P^ + 5)=0, or 



.(60) 



^ 4a/2-7(Ag-77o^) 
At ps the physical region is split into two parts i.e. matter part is confined 
m. p < ps and exterior part starts from p> Ps where a is non-negative. 

One may think that the physical singularity starts from the origin which 
makes Ps = (60). But it can not be allowed due to Coulomb's repulsive force 
of the fluid, or the fluid can only be compressed to a radius at which Coulomb 
interactions balance the the gravitational collapse effects. This restriction can 
be seen from conservation relation (15), that is time derivative of p must be 
positive. By taking the charge parameter limit highly big values (770^00) 
in /i, it gives us a relation to be satisfied by / 

Since / > and /// < for all that is / is decreasing function with 
time, then the equation is always negative or zero. Zero case corresponds 
to static solution but for the collapse situation wc get only negative values. 
Thus, existence of the electric charge violates the conservation of the energy 
momentum for p — . In other words. The conservation of the energy- 
momentum does not allow formation of the shell focusing singularity. For 
uncharged fluid 770 = 0, the central singularity can be reached and / = 
^J5Xq/A(5 gives us shell focusing singularity formation period /. 

For example, if we take time dependent function for neutral matter as 
/ = Q-i'^-V^)^ , c > with (3 — ^ — in (60) , the pressure will be singular 
at the radial distance ps ~ f"^ — e^'^~^^°^^^ * which means singularity will 
form later than Aq free case but the radius of the singularity surface will be 
greater than that of Aq. 
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Spacetime singularities. 

The Kretschmann scalar which is the square of the Riemann tensor and 
defined hy K — R^uXa R^^^'^ gives the essential, coordinate independent sin- 
gularities of the spacetime. For the spherically symmetric isotropic spacetime 
it can be given as 

Qb JjS p2 L 



K 



-16aH^pa'b{2apbb' - b'^a) + %a%^ pa'b' {2a' - pb") + 32a^b^pb'b(pb' - b) 
+ I2ab^p%{ab - 2ab) + 8a%^pb"{p + 2b') + 4a%^p^a"{aa" - 2bb) 
-16a%p%"{a%'^ + b%^) - Sa^^pba'ipb' + 26a') + 16a%%p%' {2ba' - ab') 

(62) 

K has polynomial singularities in a, 6, p and divergent as ^/gu = a — >■ 
(physical singularity), or p^O,R = bp ^0 (central singularity) . 

If the future directed non-spacelike (timelike or null) curves terminate in 
the past at the singularity then the singularity is called naked otherwise it 
is covered. The procedure is coordinate-free method. To clarify the nature 
of the singularities the future directed non-spacelike geodesies are examined, 
specially null geodesies [2, 14] . Outgoing radial null geodesies of the isotropic 
spacetime (4) are given by 

dp a 

If null geodesic equation is written in terms of physical radius R and u — p"" 
we get 

'du ~ ap^ \ ^ aR/J ' ^ ' 
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The singularity is naked if the null geodesies terminate in the past at the 
singularity with positive finite value and it is hidden or covered if the limit 

hm — = hm 1 + -— (65) 

p-^ps du p-^Psap'^-^y aR'J ^ ^ 

is negative. As stated before, the existence of the electric charge restricts the 
formation of the central singularity, only uncharged matter collapse {rio = 0) 
ends with central singularity. If all constants about the fluid defined in (24) 
are non-zero except 770 = and for a 7^ singularity starts from origin. The 
second term in the parenthesis becomes zero as — > since R' > and 
R — ph. Therefore, limit becomes positive and equal to 1, then the central 
singularity is naked. If we take r]Q — (5 — ^ — Q in. (60) at which a = 
singularity forms at Ps — ^J^X^/AaP . In this limit the (65), therefore the 
nature of the singularity becomes parameter dependent. 

In this case the limit becomes related to the sign of the time derivative of 
/ and the constant a as sgn(f)/sgn(«) . For (a > 0, / > 0) or (a < 0, / < 
with 5 < 0), limit becomes positive and therefore the central singularity 
becomes naked. But, for (a > 0, / < 0) or (a < 0, / > with 5 < 0) 
limit is negative therefore, the central singularity is covered. It is possible 
to give a lot of example such / functions so that sgn/ changes with time 
during process other than / = ^ct-^fK^t^ where "bounce" situation occurs. 
Otherwise the process will be collapse or expansion. 

Furthermore, it should be also emphasize that the physical singularity 
(a ^ singularity) coincides with the apparent horizon for the extremal 
case Aq = ?7o and /9 = . Apparent horizon is the boundary surface of the 
trapped regions and makes (51) zero 

A^b-^R'^ -a-^R^ ^0. (66) 
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5 Conclusion 



In this work spherical symmetric charged solution of Einstein field equations 
is given in the presence of cosmo logical constant. The matter is considered as 
ideal fluid and subjected to the weak energy condition. Two speciflc examples 
of this type of solution isotropic RN-(anti)de Sitter and charged Mc Vittie- 
de Sitter solutions are given. When the matter is confined in a region, the 
exterior spacetime exterior is taken RN-de Sitter and to complete analysis 
the matching conditions are examined. In these calculations, in the name 
of the generality, we started with the different cosmological constant for the 
interior Aq and exterior A regions but the junction conditions, continuity of 
the energy-momentum tensor in the radial direction gives their equality. 

In the reference [9], spherical symmetric gravitational collapse of charged 
fluid is studied in black hole formation point of view, that is, strong cosmic 
censorship hypothesis is considered. In this work, by considering weak cos- 
mic censorship hypothesis we see that besides formation of the black hole 
the process allows formation of the naked singularity and the existence of 
the cosmological constant permits bouncing situations as well. The singular- 
ity structure analysis is done by using coordinate free null geodesic method. 
The initial data about the matter (energy density, pressure and the elec- 
tric charge) therefore, constants about the solution determine the final fate 
of the collapse. Existence of the electric charge prevents formation of the 
central singularity and the cosmological constant causes bouncing situation 
for both charged and neutral matters. Uncharged matter distribution allows 
formation of the central singularity and it will be naked. 

The results are compatible with the results given in the literature. For 
further studies, the gravitational collapse phenomena can be studied for the 
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spacetimcs other than isotropic form. 
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